New predictions for inclusive heavy-quarkonium P-wave decays 
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We show that some NRQCD colour-octet matrix elements can be written in terms of 
(derivatives of) wave functions at the origin and non-perturbative universal constants once 
the factorization between the soft and ultrasoft scale is achieved by using an effective field 
theory where only ultrasoft degrees of freedom are kept as dynamical entities. This allows 
us to derive a new set of relations between inclusive heavy-quarkonium P-wave decays into 
light hadrons with different principal quantum number and with different heavy flavour. In 
particular, we can estimate the branching ratios of bottomonium P-wave states by using 
charmonium data. 



Inclusive P-wave decays to light hadrons have 
proved to be an optimal testing ground of our un- 
derstanding of heavy quarkonia. The use of NRQCD 
jl],D allowed a description of these decays in terms of 
expectation values of some 4-heavy-quark operators 
at a quantum-field level in a systematic way. Besides 
the so-called colour-singlet operators, for which their 
expectation values could be related to wave func- 
tions in an intuitive way, there were also colour-octet 
operators. The latter were decisive in solving the 
infrared sensitivity of earlier calculations ||. It has 
been thought so far that these colour-octet expecta- 
tion values could not be related with a Schrodinger- 
like formulation in any way. 

We show in this letter that it is not so. For cer- 
tain states, the expectation values of colour-octet 
operators can also be written in terms of wave func- 
tions and additional bound-state-independent non- 
perturbative parameters. We shall focus on the op- 
erators relevant to P-wave decays into light hadrons, 
but it should become apparent that this is a general 
feature. 

The line of developments that has led us to this 
result is the following. It was pointed out in Ref. 
j| that NRQCD still contains dynamical scales, 
which are not relevant to the kinematical situation 
of the lower-lying states in heavy quarkonium (en- 
ergy scales larger than the ultrasoft scale: mi) 2 , v 
being the relative velocity of the heavy quark and 
m its mass). Hence, further simplifications occur if 
we integrate them out. We call pNRQCD the re- 
sulting effective field theory (as in ||; note that in 
||, in the situation Aqcd ^ mv 2 , the EFT was 
called pNRQCD'). When the typical scale of non- 
perturbative physics, say Aqcd, is smaller than the 
soft scale mv, and larger than the ultrasoft scale 



m» 2 , the soft scale can be integrated out pertur- 
batively. This leads to an intermediate EFT that 
contains, besides the singlet, also octet fields and 
ultrasoft gluons as dynamical degrees of freedom 
j|,D. These are eventually integrated out by the 
(non-perturbative) matching to pNRQCD || . When 
Aqcd is of the order of the soft scale, the (non- 
perturbative) matching to pNRQCD has to be done 
in one single step. This framework has been devel- 
oped in a systematic way in Ref. [^) . 

In this letter we will compute the inclusive P-wave 
decay widths into light hadrons at leading order for 
Aqcd mv 2 by using pNRQCD. In this situation 
the singlet is the only dynamical field in pNRQCD 
(Goldstone bosons are also present, but they play 
a negligible role in the present analysis and will be 
ignored) , if hybrids and other degrees of freedom as- 
sociated with heavy-light meson pair threshold pro- 
duction develop a mass gap of O(Aqcd), as we will 
assume in what follows |6]j^] , or if they play a minor 
role in the heavy-quarkonium dynamics. Therefore, 
the pNRQCD Lagrangian reads |||| 

A>NRQCD = Tr|st(^-/i)s|, (1) 

where h is the pNRQCD Hamiltonian, to be deter- 
mined by matching the EFT to NRQCD. The total 
decay width of the singlet heavy-quarkonium state 
is then given by 

r = -2 Im (n, L, S, J\h\n, L, S, J), (2) 

where \n,L,S,J) are the eigenstates of the Hamilto- 
nian h. The imaginary parts are inherited from the 
4-heavy-fermion NRQCD Wilson coefficients and, 
for P-wave decays, first appear as local (delta-like) 



C(l/m 4 ) potentials in the pNRQCD La grangian. 
The relevant structure reads (we shall concentrate on 
potentials, which inherit imaginary parts from the 
NRQCD operators and which contribute to P-wave 
states at first order in quantum-mechanical pertur- 
bation theory (QMPT)): 
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where Tgj corresponds to the spin and total angular 
momentum wave-function projectors. What is now 
left is to compute F$j, i.e. to perform the matching 
between NRQCD and pNRQCD. For the situation 
A): when mv » Aqcd 3> mv 2 , by taking the results 
of ||, and for the more general situation B): when 
Aqcd < mv, by using the formalism of Refs. ||. In 
both situations we get: 



F SJ =-2N c Jmf 1 ( 2S+1 Pj) 
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where /i( 25+1 Lj) and /8( 25+1 L,/) are the short- 
distance Wilson coefhcients of NRQCD as defined 
in Ref. § and 

£=T F drr 3 <«?E Q (T,O)$ ab (r,0;O). g E b (r,O)). (5) 
Jo 

A) P-wave potentials for mv 3> Aqcd 3> mv 2 . 
In this case the matching from NRQCD to pNRQCD 
at the scale Aqcd can be done in two steps. In the 
first step, which can be done perturbatively, we in- 
tegrate out the scale mv and end up with an EFT, 
which contains singlet (S) and octet (O) fields as dy- 
namical degrees of freedom. At the next-to-leading 
order in the multipole expansion the Lagrangian 
reads |,|: 

C = Tr {s f (id -h s )S + +O f (<D - K) o} 

rr, „ n , Otr- ff EO OtOr-gE] 
+Tr <^ O f r • gE S + h.c. + ^ + — - 1 
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where h s and h have to be determined by matching 
to NRQCD. They read as follows 



V 2 

hs = -^+V s (r)- 
m 

+ N c f 1 ( 2S+1 F J )T£- 



= -— + V (r)+... 
m 

-T F f 8 ( 2S ^S s )T s S -^fi 
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neglecting centre-of-mass recoil terms; T$ corre- 
sponds to the total spin projector. Beyond 0(l/m) 
we have only displayed the terms that are relevant to 
our calculation. In the second step we integrate out 
(non-perturbatively) the gluons and the octet field 
ending up with the pNRQCD Lagrangian ([!]). The 
Hamiltonian h has to be determined by matching the 
two effective field theories. It reads h = h s + 8h s , 
with (at leading non- vanishing order in the multipole 
expansion) 



Otic = —I 
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where consistency with Aqcd 3> mv 2 requires an 
expansion of the exponentials of h D and h s . Taking 
into account that we are interested in P-wave states, 
only the perturbation that puts one h a to each side 
of the 0(l/ra 2 ) S-wave potential survives at leading 
order. The final result reads: 
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which plugged into Eq. (||) gives Eq. (Q) . This shows 
how a colour-octet operator in NRQCD becomes a 
colour-octet potential in the EFT of Eq. (Q) and, 
eventually, contributes to a colour-singlet potential 
in pNRQCD, which is one of our main points. 
B) P-wave potentials for Aqcd < mv. 
In the case Aqcd < mv the matching from NRQCD 
to pNRQCD at the scale Aqcd has to be done di- 
rectly, since no other relevant scales are supposed to 
lie between m and mv. The only dynamical degree of 
freedom of pNRQCD is the heavy-quarkonium sin- 
glet field S. The Lagrangian has been written in 
(|lj). The Hamiltonian h is obtained by matching 
(non-perturbatively) to NRQCD, order by order in 
1/m, within a Hamiltonian formalism ||. In this 
letter we only sketch the main steps of the deriva- 
tion. In short, we can formally expand the NRQCD 
Hamiltonian in 1/m: 



#NRQCD - #NRQCD + ~ #NRQCD H ■ ( 10 ) 

The eigenstates of the heavy-quark-antiquark sector 
can be labelled as 

|.g;xi,x 2 ) = |g;xi,x 2 ) (0) + — \g; x 1} x 2 ) (1) H , 

— — m — 

where g labels the colour-related degrees of freedom 
(we do not explicitly display spin labels for simplici- 
ty). Assuming a mass gap of O(Aqcd) much larger 
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than 



all the excitations (g ^ 0) decouple and 



the ground state (g — 0) corresponds to the singlet 
state. Therefore, the matching condition reads 

(0;xi,x a |fr|0;xi,34) (11) 
= /i(xi,x 2 , V X1 , V X2 )<5( 3 )( Xl -xi)5( 3 )(x 2 -x 2 ). 

Up to 0(l/m A ) the imaginary contributions are only 
carried by the Wilson coefficients of the dimension 6 
and 8 4-heavy-fermion operators in NRQCD. Since 
we are only interested in Eq. (^|) a huge simpli- 
fication occurs and only two contributions survive. 
From the dimension 8 operators we obtain 



Im5h5^\x x - x[)5^(x 2 - x' 2 ) 
= — A Im (0) (0; Xl , x 2 |f$> |o ; x ; , x ' 2 )(0) 



(12) 
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On the other hand, we also have contributions from 
the iteration of lower-order 1/m corrections to the 
NRQCD Hamiltonian with the dimension 6 4-heavy- 
fermion operators. The only term that contributes 
to Eq. (|) is 



Im5h5^(xi - x' x )5^(x 2 - x' 2 ) 



(13) 



Tm( 1 )(0;x 1 ,x 2 |ff^ QCD |0;x' 1 ,x 2 ) 



The explicit computation of the right-hand side of 
Eq. ([l3]) gives (as far as the P-wave contribution is 
concerned) Eq. <^). Therefore, the sum of the con- 
tributions from Eqs. ( |l2| ) and (O) coincides with 
Eq. (|3|), after the replacements (H) and (||). 

We can now obtain the decay widths by using Eq. 
(|). At first order in QMPT, we obtain: 
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where x l QJ {nP) ~ XQj(nP), X Q QJ (nP) := h Q (nP) 
(Q = b, c), n is the principal quantum number, 
and i?Qni(r) the radial wave function at leading 
order. Comparing with the standard NRQCD for- 
mula, where spin symmetry has already been used, 
we have: 



{h Q {nP)\O % { l S )\h Q {nP))^) = 
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T F £(n). 
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The information gained with this formula is that 
all non-perturbative flavour and principal quantum 
number dependence is encoded in the wave func- 
tion, as in the colour-singlet operators. The addi- 
tional non-perturbative parameter £ (/x) is universal: 
it only depends on the light degrees of freedom of 
QCD. This implies that the following relation be- 
tween decay widths is also universal: 
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It is interesting to notice that the UV behaviour 
of £ has the logarithmic divergence: 



£(/x) ~ 12 N C C F — In 

7T 



(17) 



which matches exactly the IR log of the 0(a s ) cor- 
rection of Im/i( 2 Pj), and hence the cancellation 
originally observed in Uj] is fulfilled. Then one could 
consider the LL RG improvement of £ by using the 
results of Ref. S for the running of the octet-matrix 
element. One obtains (/3q = 11CU/3 — 4n/Tp/3): 



£(fi) = £(//) 
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Let us apply the above results to actual quarko- 
nium, under the assumption that our framework, 
discussed in the paragraph before Eq. ([!]), provides 
a reasonable description for the P-wave states ob- 
served in nature. The numerical extraction of £ is a 
delicate task, since several of the Wilson coefficients 
(see Ref. for a full list of them) have large next- 
to-leading order contributions, which may spoil the 
convergence of the perturbative series. This problem 
is not specific of our formalism, but belongs to the 
standard formulation of NRQCD. Here, in order to 
give an estimate, we only use those data that pro- 
vide more stable results in going from the leading to 
the next-to-leading order, more precisely the aver- 
age of Eq. © for (J, S) = (1, 1), (J', S') = (0, 1), 
and (J,S) = (1,1), (J',S") = (2,1). The experi- 
mental data have been taken from |8| and updated 
accordingly to P,10) . Our final value reads 



5(lGeV)=5.3^(exp), 



(19) 



where we have used the NLO results for the Wilson 
coefficients with a LL improvement. The errors only 
refer to the experimental uncertainties on the decay 



3 



widths. Theoretical uncertainties mainly come from 
subleading operators in the power counting (0(v) 
suppressed) and subleading terms in the perturba- 
tive expansion of the Wilson coefficients (0(a s ) sup- 
pressed), whose bad convergence may affect consid- 
erably the figure of Eq. @. We feel, therefore, 
that further studies, maybe along the lines of Refs. 
fllf , are needed before a complete numerical analy- 
sis, including theoretical uncertainties, can be done. 
In any case, the above figure is compatible with the 
values that are usually assigned to the NRQCD octet 
and singlet matrix elements (e.g. from the fit of ju| 
one gets £{1 GeV) = 3.6lJ'g(exp)). The above figure 
is also compatible with the charmonium (quenched) 
lattice data of Jl3], whereas, if the running ( |l8| ) 
is taken into account, bottomonium lattice data, 
quenched Jll| and unquenched |Q, appear to give 
a lower value. Note that, in the language of Refs. 
@0, Eq. © reads £(/x) = 81ml c H s (fi)/Hi\ blC , 

which implies H^(ij) /TL\\b x Hi/Hs(^)\ c = ml/ml- 
For all quarkonium states that satisfy our assump- 
tions this equality must be fulfilled by lattice results 
for any number of light fcrmions and for any value 
of the heavy-quark masses. 

By using the estimate (|l9|), we can also predict 
the branching ratios for the n = 1,2 P-wave bot- 
tomonium states. We obtain 



r(xj x (iP)) _ rfrk(2P)) 
r(xUiP)) r(xL(2P)) 



= 0.50 



+0.06 



(20) 



where only the errors inherited from Eq. ( |l9| ) have 
been included. For what concerns theoretical uncer- 
tainties, the comments after Eq. ( |l9| ) apply also here 
(with a better behaviour of the perturbative series). 
Note that the first equality holds independently from 
Eq. ( p^ ) and from the use of charmonium data and, 
hence, provides a more robust prediction. The re- 
maining branching ratios can be obtained using spin 
symmetry. Notice also that, although no model- 
independent predictions can be made for the decay 
widths (they depend on the wave function at the 
origin, which is flavour and state dependent), our 
results allow any model that gives a definite value 
to -Rq k1 (0) to make definite predictions. 

In conclusion, we have exploited the fact that 
NRQCD still contains irrelevant degrees of freedom 
for certain heavy quarkonium states, which can be 
integrated out in order to constrain the form of the 
matrix elements of colour-octet operators. We have 
focused on the operators relevant to P-wave decays, 
which allowed us to produce concrete, new, rigorous 
results. However, it should be clear from the struc- 
ture of the pNRQCD Lagrangian itself, that similar 
results can be obtained for matrix elements of any 



colour-octet operator. 
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